arXiv:1504.04601vl [cond-mat.mtrl-sci] 17 Apr 2015 


Electrically Engineered Band Gap in 
Two-Dimensional Ge, Sn, and Pb: A First-Principles 
and Tight-Binding Approach 


Thaneshwor P. Kaloni,*’’*' Mohsen Modarresi,^ Muhammad Tahir,Mahmood 
Rezaee Roknabadi,^ Georg Schreckenbach,'*' and Michael S. Freund'*' 

Department of Chemistry, University of Manitoba, Winnipeg, MB, R3T 2N2, Canada, Department 
of Physics, Ferdowsi University of Mashhad, Mashhad, Iran, and Department of Physics, 
Concordia University, 7141 Sherbrooke Quest Montreal, Canada H4B 1R6 

E-mail: thaneshwor.kaloni(g)umanitoba.ca,-i-1 -204-952-2900 


*To whom correspondence should be addressed 

^Department of Chemistry, University of Manitoba, Winnipeg, MB, R3T 2N2, Canada 
^Department of Physics, Ferdowsi University of Mashhad, Mashhad, Iran 

“Itoepartment of Physics, Concordia University, 7141 Sherbrooke Quest Montreal, Canada H4B 1R6 


1 




Abstract 


First-principles calculations were performed to investigate the electronic structure of two- 
dimensional (2-D) Ge, Sn, and Pb without and with the presence of an external electric field 
in combination with spin-orbit coupling. Tight-binding calculations based on four orbitals per 
atom and an effective single orbital are presented to match with the results obtained from first- 
principles calculations. In particular, the electronic band structure and the band splitting are 
investigated with both models. Moreover, the simple k ■ p model is also considered in order 
to understand the band splitting in the presence of an external electric field and spin-orbif 
coupling. A large splitting is obfained, which is expecfed fo be useful for spinfronic devices. 
The fair agreemenf befween fhe firsl-principle, k ■ p model, and fighf-binding approaches leads 
fo a fable of paramefers for fulure fighf-binding sfudies on hexagonal 2-D nanosfrucfures. By 
using fhe fighf binding paramefers, fhe fransporf properties of fypical 0-D friangular quanfum 
dofs befween fwo semi-infinife elecfrodes in fhe presence of spin-orbif coupling are addressed. 
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INTRODUCTION 


After the diseovery of graphene, many two-dimensional (2-D) materials have been predieted, 
in particular, materials composed of group IV elements.-"^ These materials are very interesting 
not only for fundamental research but also for electronic device applications.- For example, H- 
passivated graphene called graphane,-^ a silicon analog of graphene called silicene,- and a germa¬ 
nium counterpart of graphene called germanene- have attracted the interest of many researchers. 
Silicene and graphene are relatively similar in nature due to the fact that Si and C belong to the 
same group IV in the periodic table. The sp^ hybridization is energetically favourable in case of 
Si, while sp^ hybridization is energetically favourable for graphene. Experimentally, silicene has 
been synthesized on metallic substrates, such as Ag,—ZrB 2 ,— and Ir(lll).— Moreover, theo¬ 
retically it has also been studied on semiconducting substrates, such as h-BN and SiC““-^ and it 
has been found that the semiconducting substrates could potentially be a way for obtaining free¬ 
standing silicene. First-principles calculations have shown that silicene is particularly interesting 
due to the formation of graphene-like dispersion.—"— 

A rather weak interaction of germanene with the semiconducting GaAs(OOOl) substrate has 
been demonstrated in such a way that a quasi-free-standing germanene could be realized. With¬ 
out H passivation of the dangling bonds of Ga terminated GaAs(OOOl), the support strongly in¬ 
teracts with the substrate, while the interaction can be strongly reduced by H intercalation at the 
interface of GaAs(OOOl) and germanene. Recently, germanene has been synthesized on a Pt(lll) 
surface.— The grown germanene exhibits a buckled structure where a 3 x 3 x 1 supercell of 
germanene coincides with a \/T9 x \/T9 x 1 supercell of the Pt(l 11) surface. The structural, elec¬ 
tronic, and magnetic properties of germanene can be influenced strongly by doping. For example, 
the Quantum Anomalous Hall effect has been predicted for V-doped germanene^ and the Quan¬ 
tum Spin Hall effect— has been predicted for H-doped germanene. Furthermore, the application 
of mechanical strain also influences greatly the structural properties of germanene.—*^ 

Since the ionic radii of the elements increase from Ge to Pb, they increasingly promote sp^ 
hybridization as a result. The magnitude of the buckling in 2-D C, Si, Ge, Sn, Pb is found to be 0.00 
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A, 0.46 A, 0.68 A, 0.84 A, and 0.89 A, respectively.— In fact, the buckling plays a crucial role in the 
engineering of the electronic structure (band gap) in these materials.— Spin-orbit coupling 
is strong in these novel 2-D materials, especially for the heavier atoms. Thus, it is important to 
incorporate it into the theoretical studies in order to achieve accurate results. Materials with strong 
spin-orbit coupling are potential candidates for applications in spintronic devices.— As already 
mentioned, these materials are buckled and they produces a very small band gap. The buckled 
structure is comparatively favorable for the application of a perpendicular electric field, which in 
fact allows tuning the band gap easily.— The band gap, in turn, is essential to operate electronic 
devices such as transistors. Recently, silicene-based transistors have been realized experimentally, 
which can be operated at room temperature.— Thus we are expecting interesting results for the 
systems under study. 

The magnitude of the spin-orbit coupling is found to be 1.10~^ meV, 4.2 meV, 11.8 meV, 36.0 
meV, and 207.3 meV, respectively for 2-D C, Si, Ge, Sn, Pb.— The combined effects of the electric 
field and spin-orbit coupling for 2-D Ge, Sn, and Pb have not been addressed so far. To addressed 
this issue, density functional theory (DFT) based calculations have been performed and the ob¬ 
tained results were verified with the help of the tight-binding (TB) approach.— In particular, the 
electronic structure (the nature of the band gap) and the band splitting were investigated. Further, 
the k- p model— was used to understand the band splitting. Excellent agreement between the DFT 
and TB calculations was found. The obtained TB parameters can therefore be used for the inves¬ 
tigation of larger nanostructures. As an example and proof of principle, we address the effects 
of spin-orbit coupling on the transport through a 0-D triangular quantum dot (TQD) between two 
semi-infinite electrodes. 
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METHODOLOGY 


DENSITY FUNCTIONAL THEORY (DFT) CALCULATIONS 

DFT calculations have been performed with the generalized gradient approximation (PBE)— as 
implemented in the Vienna Ab-initio code (VASP).— A plane wave eutoff energy of 450 eV and a 
Monkhorst-Pack k-grid of 24 x 24 x 1 have been employed.— A 1 x 1 x 1 unit eell of 2-D Ge, Sn, 
and Pb with a lattice eonstant a of 4.06 A, 4.67 A, and 4.93 A respectively, has been taken from 
Ref.^^ Moreover, a finite external eleetrie field has been applied by using the approaeh deseribed 
in Refs.'^ 

TIGHT-BINDING APPROACH 

To verify the results obtained from the DFT calculations, a TB Hamiltonian has been used for a 
unit eell that contains two atoms by following Ref.— Only the hopping between nearest neighbors 
in the unit eell was considered. The sp^ basis has been employed to expand the Hamiltonian. As a 
result the total Hamiltonian of the unit cell beeomes a 8 x 8 matrix, which is given by 

^0 = 52 + 52 ^iJ,a,l5Cla^o^i,a,cT- ( 1 ) 

i,a,a <i,j>,aP,G 

Here, ^ and Q a,c 7 represent the electron creation and annihilation operators in the hexagonal 
lattiee sites for eleetrons with spin (a) and orbital (a), respectively and <ij> represents pair of 
nearest neighbor. The first term of eq” (1) is the on-site energy for a of the i^^' atom, and the seeond 
term deseribes the hopping between nearest neighbor atomie sites < 6 j > for different values of 
a. The hopping parameters were adopted from the Slater-Koster method.— 

The Hamiltonian in the reeiproeal space is obtained by performing a Fourier transformation of 
the real-spaee Hamiltonian. The band structure is given by the k-dependent eigenstates of the trans¬ 
formed Hamiltonian matrix. Due to the buckling in the atomic structure of the hexagonal lattices, 
the presenee of a statie external eleetrie field induees a voltage differenee between the two atomie 
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positions in the unit cell. Assuming that the electric field is homogeneous along the z-direction, the 
electric field included Hamiltonian can thus be written as, = E Y.i,a,aZiCj(jCi,a,a, where Zi is 
the z-component of the atomic position of the atom Therefore, the effect of an external electric 
field is similar to a site dependent on-site energy for each atom. In principle, the buckling of the 
hexagonal lattices enhances the overlap between n and o states, and thus, stabilizes the system. 
As a result, one should consider at least four orbitals per atom to study these nanostructures. 

The interaction between electron spin and orbital angular momentum couples both degrees of 
freedom to each other. The spin-orbit coupled Hamiltonian can be expressed as. 


H 


soc 


= XL.S = l 




( 2 ) 


where A is the strength of the spin-orbit coupling, and L± and are the ladder operators for 
orbital and spin angular momentum, respectively.— In the presence of both spin-orbit coupling 
and an external electric field, the total Hamiltonian becomes = Hq-\- . By inclusion 

of the spin degrees of freedom, the total Hamiltonian becomes a 16 x 16 matrix. For the sake of 
simplicity, the overlap matrix is taken as an identity matrix and the orbital dependence of the 
strength of the spin-orbit coupling is ignored. 

An effective Kene and Mele Hamiltonian is also considered to describe the hexagonal lat¬ 
tices.—’^ It can be expressed as. 


<ij> «ij» 

Here, s, is the z-component of the Pauli matrix and and are the effective hopping 
between nearest neighbor sites and spin-orbit strength for next nearest neighbor sites, respectively, 
while = ±1, -|-/— represents clockwise/anti-clockwise orientation of the hybrid orbitals at the 
two nearest neighbor sites i and j. 

By using the effective model the electric transport is studied through a TQD between two semi¬ 
infinite structure-less electrodes as shown in Fig. 1. The TQD has localized zero-energy states on 
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Figure 1: Symbolic representation of the top view of a TQD between two semi-infinite electrodes, 
where is the number of atoms on the zigzag edge, which determines the size of TQD. 

the edges, which are split up or down in the presence of a Hubbard term by inducing magnetic 
moments on the zigzag edges.— The GreenaAZs function approach— is a well-known method 
for calculating the transport properties through a nano-structure between two leads. Here, a TQD 
between electrodes has been studied in the presence of intrinsic spin-orbit coupling. 

In order to describe the above system, one should consider a Hamiltonian that contains the 
TQD, the left and right electrodes and the coupling between them. The net effect of the two 
semi-infinite electrodes is described by using self-energy terms that can be calculated by a simple 
iterative method. Here, for simplicity the wide-band approximation is used which ignores the 
real part of the self-energy and estimates the imaginary part with a constant value for all energy 
intervals.— The self-energy operator is useful in calculating the retarded GreenaAZs function, 
defined as, G''(e) = {{£ -\- id^)l — Hjqd — — . In this equation, El//? the self-energy 

operators for the left and right electrodes, I is the identity matrix and Htqd is the Hamiltonian 
of the quantum dot. The transmission spectrum as a function of energy through the quantum dot 
is calculated by using the Fisher-Lee relation— r(e) = Trace [r/^G'‘(e)r^G'‘(e)], where = 
is the conductance of the noninteracting electrons. The current-voltage curve is calculated 
by using the Landauer-Buttiker equation and the transmission spectrum.— 
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Figure 2: The eleetronie strueture of 2-D Ge for different values of the external eleetrie field E^. 
Only the zoomed area around the K-point is shown for (a-d). The band splitting is represented by 
a pink square. 


RESULTS AND DISCUSSION 

ELECTRONIC STRUCTURE 

Table 1: The magnitude of the external eleetrie field (in V/A), band gap Egap (in meV), and band 
splitting (in meV) for 2-D Ge, Sn, and Pb, where m represents metallie states and * represents 
an indireet band gap. 


System 

Ge 

Sn 

Pb 

Ez 

F 

^gap 

Ap 

E 

gap 


F 

^gap 

A 

sp 


DFT 

TB 

DFT 

TB 

DFT 

TB 

DFT 

TB 

DFT 

TB 

DFT 

TB 

0.0 

22.9 

34.2 

0.0 

0.0 

73.4 

77.4 

0.0 

0.0 

409.1 

366.3 

0.0 

0.0 

0.001 

20.1 

33.7 

3.5 

0.45 

76.3 

76.8 

0.71 

0.6 

408.5 

365.8 

0.4 

0.5 

0.005 

23.8 

31.5 

23.1 

2.7 

78.2 

75.5 

0.76 

2.9 

406.1 

363.8 

4.5 

2.5 

0.01 

47.4 

58.7 

24.5 

5.4 

79.5 

71.3 

7.51 

6.1 

404.9 

361.3 

5.4 

5.0 

0.5 

31.5 

242.2 

25.8 

35.2 

309.9 

227.8 

76.6 

77.9 

157.2 

117.1 

266.9 

255.2 

1.0 

115.4" 

518.7" 

26.1 

33.6 

318.3 

530.1 

81.2 

79.2 

m 

124.2 

0.0 

359.2 

1.5 

408.8" 

792.2" 

34.8 

37.7 

328.7 

826.9 

86.2 

81.3 

m 

352.4 

0.0 

357.0 

2.0 

570.4" 

1063" 

43.3 

39.3 

228.3" 

1117.2" 

90.1 

84.1 

m 

564.0 

0.0 

354.3 


The eleetronie strueture of 2-D Ge, Sn, and Pb was studied without and with the presenee of an 
external eleetrie field by employing DFT ealeulations and the TB approaehes. Spin-orbit eoupling 
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has been taken into account in all the calculations. The calculated band structure of 2-D Ge is 


addressed in Fig. 2. The calculated band gap (E^ap) in the absence of an electric field is found to 
be 22.9 meV with no band splitting, which is in good agreement with the available reports.— 

The values of and Asp increase with increasing magnitude of the electric field.— The A^p is 
increased because of the combined effect of spin-orbit coupling and an external electric field, which 
agrees well with a previous report for silicene.— Whereas, the band splitting can be attributed to 
the fact that the external perpendicular electric field removes the space inversion symmetry, the 
absence of space inversion symmetry lifts the spin degeneracy due to spin-orbit coupling; as a result 
the splitting between states has been observed.— It is worth mentioning that transition from a direct 
band gap (at the K-point) to an indirect band gap (the valence band maximum and conduction band 
minimum at the K- and E-points of the Brillouin zone, respectively) is obtained above an external 
electric field of 1.0 V/A. A direct to indirect band gap transition occurs under an external electric 
field because the electric field enhanced spontaneous polarization. This polarization is responsible 
for the direct to indirect band gap transition.In addition, a perpendicular electric field reduces 
the band gap and thus, the transition is possible. The data sets are summarized in Table 1. 

Eor low values of the electric field the electronic gap is related to the spin-orbit coupling but for 
higher values of the external field the gap results from the difference of the electric potential on 
the different atomic sites. An unscreened uniform electric field has been considered; as a result, 
for high values of the electric field strength, the tight binding model predicts a relatively bigger 
electronic gap with respect to the DET results. 

The transition from a direct to an indirect Egap has been observed previously either via strain 
or via the application of an electric field in 2-D transition-metal dichalcogenides, such as M 0 S 2 
and In addition, a large value of about 43.3 meV is obtained in 2-D Ge for an 

electric field of 2.0 V/A, which is expected to be useful for spintronic devices.-^^*^ In general, due 
to the buckled structure of 2-D Ge, the presence of an external electric field induces a difference 
in the electric potential on the two atoms of the unit cell. The potential opens an energy gap at the 
K-point which increases with increasing magnitude of the electric field. The spin-orbit coupling 
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together with the eleetrie field splits the up and down states. 


(a) = 0.0 eV/A (b) = 0.001 eV/A (c) E^ = 0.005 eV/A 
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Figure 3: The eleetronie strueture of 2-D Sn for different values of the external eleetrie field 


The ealeulated eleetronie band strueture of 2-D Sn is shown in Fig. 3. A Egap of 73.4 meV 
with Ajp = 0 is obtained in the absenee of an external eleetrie field, see Fig. 3(a). It has been 
reported that 2-D Sn also forms a Dirae-like eone without inelusion of the spin-orbit eoupling 
in the ealeulations.- The authors of Ref.^ observed a band gap of about 100 meV by removing 
the degeneraey between the bands at the F-point when the spin-orbit eoupling is switehed on. 
Therefore, it is expeeted that 2-D Sn ean host a Quantum anomalous Hall effeet,— similar to 
graphene,— silieene,— and germanene.— The magnitudes oiEgap and beeome 76.3-328.7 
meV and 0.71-86.3 meV, respeetively, for external eleetrie fields of 0.001-1.5 V/A, see Table 1. 
The transition from a direet to an indireet band gap is obtained for an external eleetrie field of 2.0 
V/A. In this ease, an indireet Egap of 228.6 meV along the F-K direetion and a A^p of 90 meV at 
the K-point and 21 meV at the F-point are obtained. The external eleetrie field results in loealized 
eigenstates and the eleetronie bands at the K-point beeome flatter. 

In Fig. 4, the eleetronie band strueture of 2-D Pb is addressed. Upon inelusion of spin-orbit 
eoupling, an Egap of 409.1 meV is obtained in the absenee of an external eleetrie field. In this ease 
Egap is deereased with inereasing field strength and the material beeomes metallie for an eleetrie 
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Figure 4: The electronic structure of 2-D Pb for different values of the external electric field 

field of 1.5 V/A or higher. It is worth to mention at this point that a semiconductor-to-metal 
transition has been demonstrated also for few layers and bulk M 0 S 2 .— An Eg^p value of 157.2 
meV with a A^p of 233.1 meV is obtained for an external electric field of 0.5 V/A, see Table 1. 
The structural parameters agree well with a previous report.— We note that a clear trend is found 
for the group IV elements, especially in the nature of the band gap. A band gap of 0.0 eV, 2.0 
meV, 22.9 meV, 73.4 meV, and 409.1 meV is obtained for C, Si, Ge, Sn, and Pb-based systems. 
The increase of the band gap with increasing atomic size can be attributed to the increase in the 
strength of spin-orbit coupling. 

In addition, to understand the mechanism of Eg^p and A^p, a two-band k- p model has been 
applied. This model is sufficient to describe the charge carriers in these novel 2-D materials with 
large intrinsic spin-orbit coupling.— An effective Dirac-like Hamiltonian is used in the vy-plane to 
describe the physics of the valence and conduction bands in the K and k' valleys as 

= v{r\d^px + dyPy) + sr\Xdz + XyOz- (4) 

Here, t] = ±1 represents the K and K' valleys, respectively, (oy, dy, d^) is the vector of Pauli 
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matrices (applicable to both the valence and conduction bands), A is the spin-orbit coupling energy, 
5 = ±1 represents the up and down spins, respectively, and v denotes the Fermi velocity of the Dirac 
fermions. Xy = llE^ is the term induced by the perpendicular electric field. It breaks the inversion 
symmetry. Finally, 21 is the perpendicular length along the buckled axis.— 

After diagonalization of the Hamiltonian of Eq. (4), the eigenvalues are obtained as 

{vhkY -\- {Xy — srjXy (5) 

Here, t = ±1 represents the conduction and valence bands, respectively. Two aspects are worth 
noticing; (1) in the limit of zero electric field, a spin-orbit coupling induced gap (Asp = X) can 
be obtained at wave vector k = 0, and (2) by increasing the perpendicular electric field, both the 
valence and conduction bands show splitting due to a combination of spin-orbit coupling and 
electric field energy (Egap = Xy — srjX), which for example is 22.9 meV, 73.4 meV, 409.1 meV 
for 2-D Ge, Sn, and Pb, respectively. With increasing perpendicular electric field, the spin-orbit 
splitting is enhanced and well resolved in both the bands. Moreover, all the results obtained using 
first-principles calculations are matched by the results from the TB approach for the lower range 
of the electric field, see Fig. 5 and Table 1. However, the TB fails to predict the properties of 2-D 
Pb in the presence of high external electric fields. 


Ge Sn Pb 



Figure 5: The band structure obtained from the TB calculations (blue lines) and effective TB 
approach (red lines) for 2-D Ge, Sn, and Pb lattices. 
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Table 2: The TB parameters for 2-D Ge, Sn, and Pb lattiees, where Eg and Ep are the on-site energy 
for s and p orbitals and all the parameters are measured in eV exeept whieh is measured in 
meV. 


System 

Es 

Ep 

Es(7 

tspo 

^ppa 

tppTt 

A 

ten 


Ge 

-5.65 

-0.8 

-1.6 

2.4 

2.1 

-1.1 

0.15 

-1.1 

2.5 

Sn 

-5.0 

0.65 

-1.2 

1.9 

1.7 

-0.8 

0.2 

-0.8 

7.0 

Pb 

-7.15 

0.8 

-0.85 

1.4 

1.4 

-0.7 

0.45 

-0.7 

35.0 


The effective TB model describes the band structure at the K-point effectively but it fails around 
the F-point. In order to fit the strength of the spin-orbit interaction, the size of the electronic gap 
and band splitting at the K-point are considered. In Table 2, the TB parameters for 2-D Ge, Sn, 
and Pb are summarized. The TB parameters are obtained according to the method of Ref.— Ea 
is the on-site energy and are different components of the hopping between nearest neighbor 
sites. The obtained A is proportional to the atomic mass of each element. By increasing the atomic 
radius the effective hopping and spin-orbit strength decreases and increases, respectively. The TB 
parameters obtained in this study could be useful for future calculations of nanostructures as will 
be illustrated in the following section. 

TRANSPORT PROPERTIES 



E-Ep (eV) 

Figure 6: The transmission spectrum of TQD between two electrodes for = 7, 15, and 21 atoms. 
In the following the transport properties are studied for TQD by using the effective tight binding 
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model. Triangular shaped graphene quantum dots have been synthesized experimentally— and it 
has been suggested that quantum dots made from the other group IV element should be triangular 
in shape as well. The ratio varies from zero for graphene up to 0.05 for Pb. For the 

self-energy, only the hopping between eleetrode and the edge atoms of the TQD is eonsidered and 
the imaginary part of the self-energy is set to 0.05 eV. The TQD is eharaeterized by the number of 
atoms on eaeh edge Nq. The transmission speetrum of different size quantum dots is addressed in 
Fig. 6. 

The transmission is symmetric and peaks are corresponding to eigenstates of the Hamiltonian. 
By increasing the size of the quantum dot, the length of the dot increases, and the transmission 
through the nanostructure decreases. In the absence of spin-orbit coupling the triangular dot has 
(Nfl — 3)/2 degenerate and localized wave functions in the zero state energy that do not contribute 
to the transmission. The spin-orbit coupling removes the degeneracy and results in localization of 
the wave functions; the transmission gap is filled with peaks. For larger dots the energy difference 
between states is decreased and the peaks in the gap become flat. The first two peaks around — 1 eV 
in the transmission curve are corresponding to the first electronic states. By increasing the strength 
of spin-orbit coupling the first transmission peaks are decreased. The applied voltage difference 
between left and right leads to an opening of the energy window and current flows in the quantum 
dot. The transmission peaks in the gap lead to steps for low voltages in the 1— V curves. By 
applying an external electric field perpendicular to the dot surface the energies of the localized 
states are shifted to lower or higher energies depending on the direction of the electric field. As a 
result the transmission spectrum becomes asymmetric with respect to the Fermi level. 

CONCLUSION 

In summary, three different 2-D hexagonal lattices, specifically Ge, Sn, and Pb, were investigated 
using first-principles (DFT) and tight-binding approaches. The electronic structure of 2-D Ge, 
Sn, and Pb without and with the presence of a perpendicular electric field in combination with 


14 


spin-orbit coupling have been addressed. The obtained results are eompared with the tight-binding 
approaeh. The nature of the band gap, a transition from semieonduetor to metal, a transition from 
direet to indireet band gap, and band splitting are found and discussed. A large splitting and tuning 
of the band gap are obtained by varying the perpendicular external eleetric field, whieh is expeeted 
to be useful for spintronic deviee applieations. We emphasize that silieene-based transistors have 
been realized experimentally already, and they can be operated at room temperature.— Thus we 
are also expeeting interesting results for our systems under study. In addition, the k ■ p model 
is eonsidered in order to understand the band splitting. Excellent agreement between the first- 
prineiples and tight-binding approaehes is found. Furthermore, the transport properties of typieal 
0-D triangular quantum dots between two eleetrodes are addressed, whieh will pave the way to the 
experimental realization.— 
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